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Honp'-rurr.-  trie  Tests  of  Independence  for  Consort’d  Data 
with  pplications  fo  H<  art  Transplant  Studi  -s 

X  X-*  XXX 

Byron  Wm.  Brown.  Jr.,  Myles  Hollander,  and  Raim  -sh  M.  Korvnr 

Abstract .  A  patient  officially  selected  as  a  heart  transplant 
candidate, in  the  Stanford  Heart  Transplant  Program,  will  receive  a 
transplant  if  he  survives  until  a  donor  is  found.  H<  will  not  reoeiv 
•  w  h*"irt  if  h  di<.'  before  a  donor  is  found,  or  if  lv  is  "der.  - 
-■  ct--d"  (’a i  bowing  notev/orthy  improvement).  For  a  candidate  who  no’ 
.ot  c*  iv-  transplant,  th  availnbl<  data  include  his  survival  t  jm» 
•  y  xu  th‘  date  i'  acceptance  into  the  program)  and,  if  d<  .  - 

.  hi  j  ■  of  ib  ;  •!<  c  ion.  For  a  candidate  who  do-  .  rec<  i  v«-  a 
i  i.narii  ,  'h  •  va  i.l  abe  un i.a  include  th>>  waiting  t  ime  for  t  donor 
te  a  "/ival  i  im  from  tie  date  of  operation.  V.  .  I  tii.g  tiim  nu 
-rvival  tim*  v  t iabF  .  may  bt  censored  at  the  closing  dati  ot  t  h< 
t  . iy .  Turnbul  1  at  rr own  111]  aim  Turnbull,  Brown,  and  Hu  [15]  con- 
ri  :cr  1  *h  a-  ‘ion:  Do-  s  cardiac  transplantation  at.  Stanford 
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prolong  life?  [ Th  y  found  there  wan  not  sufficient  evidence  to  sup¬ 
port  a  "Y  a"  answer;  in  fact  their  unulysos  showed  tint  the  posl - 
transplant  survival  times  of  transplant  patients  were  consistent  with 
pre-transplant  experience  of  all  heart  candidates.]  Here  we  investi¬ 
gate  which  variables  nr>‘  correlated  with  post-transplant  survival 
length.  Variable,  considered  are  sex,  age,  waiting  time  for  donor, 
and  date  of  transplant.  These  investigations  lead  to  new  nonpa ra¬ 
nt  trie  tests  of  independent  for  censored  data. 

1.  Introduction.  H*  art  transplantation  has  held  on  extraor¬ 
dinary  fascination  for  both  layman  and  medical  scientist  ever  since 
the  first  operation  by  Barnard  in  1967*  After  the  first  flurry  of 
operations,  two  program,  in  the  United  States  continu  d  to  accumulat' 
experience  in  tran  plantation,  one  at  Baylor  University  and  one  at 
Stanford  Univ<  rsity.  ’The  survival  experience  of  transplant  patient  : 
at  Baylor  vac  summarized  by  Mm  amor,  et  al  [  133  •  They  concluded  that 
survival  tim  of  their  transplant  pati  nts  was  ext  aided  but  the  addi¬ 
tional  life  tim*  was  not  worth  th<  attending  cost,  including  trauma; 
thf>  program  was  discontinued.  The  Stanford  xperienc<>  war,  summarized 
by  Clark,  ct  al  [2],  They  also  concluded  that  transplantation  ex¬ 
tended  life  and  the  program  at  Stanford  continues,  with  one  to  thr 
new  transplants  a  month. 

Gail  [7]  examined  th-  M'-saim  r,  et  al  [13]  and  Cl:  rk,  et  al  [2] 
papers  and  found  that  both  pipers  employed  a  n  shod  of  analysis  ‘  Kat¬ 
ie  highly  biased.  Briefly,  thf  analyses  compared  survival  tiim  s.  for 


h^art  transplant  patients  with  the  survival  times  of  transplant  candi¬ 
dates  who  died  before  a  donor  heart  was  found. 

Turnbull  and  Brown  [lU J  re-analyzed  the  Stanford  datu,  taking 
Gail's  criticism  into  account,  and,  recently,  Turnbull,  Brown  and 
Hu  [15]  hav  analyzed  the  Stanford  data  again,  summarizing  the  data  as 
of  March  1,  1973*  They  concluded  that  the  survival  times  of  heart 
transplant  patients  are  quite  consistent  with  pre-transplant  experi¬ 
ence  of  all  pre-transplant  heart  candidates  selected,  but  that  no 
d-finite  conclusion  can  be  reached  on  the  question  of  life-extending 
effects  of  transplantation,  because  there  is  too  little  long-time  pre¬ 
transplant  data  available.  However,  there  is  good  evidence  that  trans¬ 
plantation  does  not  shorten  expected  survival  time  to  any  great  ex¬ 
tent,  and  the  quality  of  life  for  some  survivors  is  remarkably 
enhanced. 

The  analyses  of  Turnbull,  et  al([ll+],  [  15 ] )  did  not  consider 
which  variables  were  correlated  with  the  length  of  post-transplant 
survival.  Here  we  address  this  question  of  correlation  considering, 
in  particular,  the  variables  sex,  age  at  transplant  date,  date  of 
transplant,  and  waiting  time  for  a  donor.  The  investigation  of  a 
possible  correlation  between  date  of  transplant  and  post-transplant 
survival  is  of  special  interest  because  it  can  be  viewed  as  an  exami¬ 
nation  of  a  trend  in  the  expected  survival  tim'1  of  transplant  patients 
(with  an  increasing  trend  possibly  reflecting  improvement  in  surgical 
techniques,  in  the  methods  of  selecting  transplant  candidates,  and  in 
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monitoring  and  caring  for  the  post-transplant  patients). 

Our  approach  is  via  Kendall's  [l?]  rank  correlation  statistic. 

This  statistic,  however,  has  to  b  modified  to  accommoda t<  right  en¬ 
coring  of  th<  dependent  variable,  typically  found  in  clinical  survival 
■*u!iec.  Two  conditionally  distribution-free  tests  ar  •  presented;  one 
(S'  ction  3)  is  analogous  to  Cohan's  [8]  adaptation  of  the  Wilcoxon  two 
sample  statistic  to  censored  data  and  the  other  (Section  4)  ic  analo¬ 
gous  to  Efron's  [6]  adaptation  of  the  Wilcoxon  two  sample  statistic  to 
censored  data.  Th-  conditionally  distribution-free  tests  of  Section;  3 
and  4  are  based  on  an  independence  assumption  concerning  the  basic 
variables  (whon  independence  is  being  tested)  and  th  censoring  vari- 
bl-  .  Thi~  assumption  [Assumption  (A) — see  Section  3]  will  not  be 
tisfi  ir.  certain  situations  of  interest.  Hence  ,  in  Section  5,  w' 
relax  this  assumption  and  provide  u  "pseudo"  conditional  test.  This 
test,  however,  do'  c  not  possess  the  distribution- fre o  property. 

Section  6  summarizes  the  applications  of  these  techniques  to  questions 
of  interest  in  th  Stanford  Heart  Transplant  Program.  Pobustn  ;ss 
comm  nts  are  raau  in  Section  r( . 

2, _ Data  d<  scription.  The  starting  date  for  a  patient  will  be 

th  date  on  which  he  is  declared  a  transplant  candidate  in  a  team 
conference,  initiating  the  search  for1  a  donor.  The  data  for  tin.  8 
candidates  ndmilt  d  to  th  program,  as  of  March  1,  19Y3?  arc  giv  >n  in 
Table  1.  For  patients  who  died  before  a  donor  was  found  for  them  and 
for  those  still  awaiting  a  heart,  we  give  their  birth  date,  date  of 
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rl  clared  acceptance  into  the  program,  sox,  and  the  number  of  day.  o-' 
survival  (to  death  or  to  the  closing  uatr  for  analysis,  March  1.  19",  ). 

■Iso  indicate  whether  the  pati"nt  is  alive  or  dead  it  the  closing 
date.  For  patients  receiving  a  new  h -art  we  give  th>  birth  date,  dot 
of  acc-  p*anc',  sex,  days  to  transplant,  days  from  transplant  to  d  >ath 
or  closing  dati  ,  and  state  (d  ad  or  aliv  )  at  closing  dal 


Tab  It  1:  Acceptance  dates,  birth  dates,  sex,  survival  times,  and 
times  to  transplant  for  Stanford  Program  patients. 
(Closing  date;  March  1,  1973) 
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9 
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12 

d 
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72 

15 
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6 

27 
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6 

72 
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12 
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. 

T 

2 
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to  tranr, 

plant 

• 

T 
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=  clay. 

from  tranr pin 

nt 

to  d 
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or  cloning  d 

S/d  * ; 

.elect,!  1 

8/21/69 

• 

S/  Du 

;el«  cteci 

2/02/7 

(a  In 

0 . 

l’or 

thin 

p  ti 

•nt  only. 

nurvival  evp,  rienc 

-  in 

kn 

own  only 

up  to 

1/1/73 

r.o 

that  t.h 

•'  T1  = 

127  < 

ntry  in 

men 

.  tired 

from 

11/1/71  is 

1/1/73 

ra  th1 

r 

th.  1  n 

3/1/ 

.  ;)■ 

J 

D  c 

1"C  'I 

/O'./  /  . 

. 

Source:  .  lor.,:  ('  .  ) 
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3.  Sinrpl  ■  scores  modification  of  Kendall's  rank  correlation 
statistic.  For  data  convicting  of  n  pairc  (X^,Y^) , . . . , (X^,Y  )  th>' 
hypothesis  of  independence  of  X  and  Y  can  be  ter.ted  ucing 
Kendall’ c  [12]  rank  correlation  statistic 

n  n 


(1) 


"  i jbi  i  ’ 
i=l  j=l  1J 


wh  re  ai>t  =  1  if  X.  >  X  ,  0  if  X±  =  Xj,  -1  if  X.  <  X.,  b  =  1 

if  Y.  ■>  Y, .  0  if  Y.  =  Y.,  -1  if  Y  <  Y..  If,  however,  either  X 

i  j '  i  j  i  j 

or  Y  (or  both)  is  censored,  we  may  be  unable  to  compute  certain  of 


fh< 


or  b’r-.  In  our  heart  transplant  studi>  c  w  >  ncounterod 


situations  vh  r  on  of  the  variables  (cay  Y)  was  right -censored; 
for  convenience  and  simplicity  we  dev-  lop  our  tests  for  this  cn.  . 
Modifications  in  situations  where  X  and  Y  ar  both  right-arid  1  ft  - 
c  nsored  arc-  notutionally  more  complex  but  can  b  d  v  loped  in  i  simi¬ 


lar  ms  nr.  r. 


When  th*  Y  variable  is  right -censored,  our  observed  paired  data 
may  b^  characterized  by  the  vector  W  =  ( (X^,Z^,5- )  •  •  •  •  *  (X^  ,7.^ ,6^  )} 
where,  for  i  =  l,...,n, 


(2)  Z^  =  minimum  (Y^,B^)  , 

and 


{  1  if 

z.  =  Y.  (that  is. 

Y. 

is  uncensored) 

( 

1  1 

i 

0  if 

-  B.  (tha4  ic, 

J-  -L 

Y. 

i 

is  c  ns, os 'd  a  t 

Known  value. 
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For  ex'  rnpl'  .  ir.  our  study  of  a  possible  correlation  between  sex 
and  length  o:  port -transplant  survival,  =  1  if  patient  i  is 
fjile  and  0  if  fcrr-le,  and  Y_.  is  the  number  of  days  patient  i  stir 

X 

viv  d  aft  r  trr  toplant.  If  patient  i  dies  before  the  closing  dat«.  . 
Y.  is  ur:c  .nsored  and  Z.  =  Y.  .  If,  however,  pati  -nt  i  is  still 

i  ii 

alive  on  th<  clo  ing  dut< ,  th. -n  Y.  is  (for  our  present  purposes) 
c  nsorecl  and  v.  in  at  d  oK  rv  Z  =  B^,  the  number  of  days  from  Lh 

Jate  of  transplant  to  the  closing  (bite.  [In  thr  notation  of  Table  1, 

■  ^  =  1  corresponds  to  th<  letter  "d,"  &  =  0  corresponds  to  "a."] 
For  th  cue-  wh'-r  Y  is  right-censored,  we  adjust  th  d-  fini- 
*  ion  of  the  b's  as  follows: 


1  if  Y.  >  Y.  , 

i  J 

0  if  Y.  =  Y.  or  if  "uncertain"  , 
i  J 

d 

-1  if  Y.  <  Y.  . 

t  J 


In  (*) ,  Y  >  Y.  (road  "Y.  is  d  finitely  greater  than  Y.")  m>. -a  ns 

i  <J  1  •) 

that,  on  the  bn  i  of  Z.,Z.,b.,&.,  we  can  infer  that  Y.  ">  Y..  Tile 

i  .J  i  J  i  J 

d 

notation  Y.  <  Y.  is  d  fined  similarly,  and  b.  .  =0  if  we  can  infer 
J  10 


th  t  Y.,  =  Y,  or  if,  because  of  censoring,  we  cannot  be  sure  of  th< 

d 

.Y.  or  • ring.  W>  have  that  (i)  Y.  >  Y.  if  (  )  .  =  1  and 

—  <j  i  J  ■  J 


1 1 

A 

or  (b)  b.  =  0, 

a  .  =  1,  and 

Z.  =  Z,.  Yi<  ha V1  (ii) 

d 

•j 

i  J 

Y. 

<  Y. 

i  (c)  e,  =  1  ■ 

ind  Z.  <  Z,, 

or  (a)  5.  =  1,  .  =  0,  ana 

X 

J 

X 

.r 

i  .) 

<y 

t~i  j 

X 

_  7 

For  (i),  g..  = 

1;  for  (ii), 

b.  ,  =  -1.  If  both  (i)  •  rut 
ij 

Th*  significance  or  the  statistic  S,  for  testing  independence , 


atm  be  or >t  •  in*  i  by  ijv<-ing  where  tin’  observed  vitlue  of  5  (G{w)  my) 

•:11c  In  ’In  }>  r;  an  Lion  distribution  obtained  by  computing  n!  v:  1 
o  om  value  for  ouch  of  the  n!  possible  ponuututions  of  tlie 

{(  ^,&^) . (.’  a  )]  obm ervntions  (ur  ping  tin-  X's  fixed).  Con- 

siuer  bli  n  th  'roup  G  of  u!  Irons  romiations  with  typical  member 


.  (V.)  =  {(Xn,.t.  ) . (X  ,Z.  .6.  ))  . 

i  '  '  1  ■  i-  i  '  v  n  l  i  ' 1  ■ 

n  11  n  n 


vdiere  (i. . i  )  is  a  permutation  of  (l,...,n).  As  a  reference 

1  n 

point  for  di  '  eimiiiin;;  the  significance  of  S(w) ,  we  consider  the 
1.  •  of  tlie  conditional  measurer; 


Or)  f(>v  =  g.  ......  (W))  «  (nO-  . 

1  11 

for  end;  pi -mutation  (i, . i  'i  o  (l . n) . 

Tla  us.  of  (^)  cun  b  Justified  i  '  tin  ind  pend  nee  of  X  nr  i 
implies  tlp/'ir.'  -pend-nci  0  X  and  (Z,5),  so  tint  when  tin  hypoth¬ 
esis  of  X.Y  i.vb  p  :  u  nc  is  true-,  or  ouch  g  c  G  ,  g(W)  would 
lwve  th'1  Gair."  distribution  an  W.  In  developing  the  conditional 
distribution -free  tests  of  Sections  3  and  4,  we  thus  take  *,*  a  basic: 
assumption 


/.  v  When  X  and  Y  nrr  irxd<  pendent, 

'  X  •  ml  (/. .  )  •  r>  imh  pendent.  . 

Whether  or  not  ^mptior.  (A)  is  .aitir.fjfd  in  practiee  m-pends  on  th 
joint  distribution  of  (X.Y.B).  [Ar,  assumption  that  implies  (a).  t>m 
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ic  unnecessarily  restrictive  is  (A'):  Wien  X  and  Y  are  independ¬ 
ent,  X,  Y.  and  B  ar*  mutually  independent.]  A  rough  (but  not 
uniformly  valid)  rule  is  that  Assumption  (A)  is  a  reasonable  assump¬ 
tion  to  make  vh  n  the  ind<  p'  ndence  of  X  and  Y  irnpli <•  e.  the  inde¬ 
pendence  of  X  and  B.  For  example,  with  =  1  if  patient  i  is 

male  and  0  if  femal  .  and  Y^  numbt  r  of  post-transplant  survival 
days  for  patient  i,  Assumption  (A)  ic  reasonable.  However,  if  we 
d<  fin  X^  to  be  the  date  of  transplant  (and  ke  p  Y  as  aboV'  ). 

then  ,,inc  B.  =  number  of  days  from  X^  to  closing  sate,  it  is  cl-  ar 

that  X  and  B  have  a  correlation  of  -1  whether  X  and  Y  ar 
independent  or  not.  (W<  can  write  B^  =  T-X^  vher  T  is  a  fixed 
constant.)  In  this  case,  even  if  X  and  Y  Sri-  independent.  X  and 

(Z,b)  ar<  dependent  and  (‘>)  is  not  justified. 

When  Assumption  (A)  is  satisfied,  our  Q- level  test  is  defii  d  as 
lows.  Let  f^^(w)  <  •••  <  S^n’^(w)  denote  the  nl  ordered 
values  o:  S(g(  ))  for  g  e  G.  Let  rp(w)  denoti  th«  probability  of 

rejecting  the  hypothesis  of  independ<  rice  when  V/  =  w.  We  set 


'  i 

if 

S(w) 

>  G(m)<vO  . 

(6) 

c?(w)  =  r(w) 

if 

e(w) 

-  s<n>(„)  , 

0 

if 

G(w) 

<  s(">(»)  , 

where  m  =  n!  -  [nl  0 ,  [x]  is  the  greatest  integer  loss  than  or 
equal  to  x  and  r(v)  is  selected  to  give  the  test  size  a. 
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Under  Assumption  (A),  the  test  defined  by  (6)  is  a  one-sided  condition¬ 
ally  distribution-free  test  of  independence  versus  the  alternative  of 
positive  association  between  X  and  Y.  One-sided  tests  versus  nega¬ 
tive  association  and  two-sided  tests  are  similarly  defined. 

For  large  n,  the  following  approximation  to  the  conditional 
distribution  of  S  can  be  used.  Under  suitable  regularity  (a  suf¬ 
ficient  condition  is  that  la., a..,  and  I.b,  ,b . , ,  each  be  of  the 

ij  ij’  ij  ij' 

order  of  n  )  the  statistic  S  is  conditionally  asymptotically  normal 
wit-h  conditional  mean 

(7)  E(S)  =  0 


and  conditional  variance 


(8) 


Var(o)  =  b[n(n-l )(n-2) J  1(Xa.  a  ,-ia?  )(fb  b  ,-Ib  ) 

i.3  i.1  i.3  i.)  ij  ij 

+  ?[n(n-l)r1(Xa?,)(r.b?,)  . 

i  J  ij 


In  the  above  summations,  each  subscript  is  summed  from  1  to  n.  The 
conditional  mean,  variance,  and  asymptotic  normality  fellow  directly 
from  Daniels  [  t  . 

In  the  special  case  when  the  distribution  of  (X,Y)  in  continuous 
and  there  is  r.o  censoring,  exact  critical  values  of  S/2,  for 
n  =  Ml)!*0,  'an  be  obtainei  from  Table  I  of  Kaarsemaker  and  von 
Wijngaarden  [ 1 0 ] ,  which  also  appear,  as  Table  A. 21  of  Hollander  and 


Wolfe  [9 ]. 

If  ve  set  X.  -  1  if  i th  patient  is  male  and  0  if  female 
then  Assumption  (A)  is  reasonable.  Vh i •  dichotomy  puts  the  problem  in 


- r. o  framework  of  the  *,wo  sample  location  problem  for  censored  data, 
for  ve  are  then  comparing  the  Y  variable  in  the  nnle  and  female 
populations .  In  ouch  a  case,  the  conditionally  distribution- free 
tent  given  by  (6)  reduces  to  Mian's  f 8 ]  procedure.  Of  course, 
Assumption  (A)  is  reasonable  for  other  types  of  X  variables  (e.g, 
physical  status  at  date  of  transplant,  or  degree  of  tissue  match 
with  donor). 

X  u  1  i  r  scores  mod  1  l'i c~ 1 1  i on  of  K-  ncbill'.:  rank  cor r •  la 

-  ar,  r  ’ .  ~  *  i .  •-  i  c  .  Th-  b’s  cl  fined  by  (a)  ignon  c-rtain  -.vailabl- 

inforawtion.  For  example,  suppose  for  the  moment,  v-  identify 

th--  post -transplant  survival  variable  T^,  of  Tab]'-  1.  Connie-  r 

patients  15  and  51.  with  7,  valuer.  1379  and  hj  reap- -ct  j v- ly .  Th- 

approach  of  Section  « cores  b. ,  =  0  sine--  b.,  -  >  ,  =  0.  Tli 

•t-Oj.e  dU 

in,  sine-  both  p-ifi  nts  were  ilive  at  th  closing  dal  -  ,  w-  nr-  un¬ 
certain  an  to  which  on-  will  accumulate  more  post -t  ran  , plant  1  i  f-  . 
However,  inc-  pnti-nt  15  La;.  b< •* -n  best  ,  to  dot-  ,  in  t.cnns  of  post- 
transplant  survival  ir.  th  Stanford  lYogram,  on-  feels  it  is  likely 
that  Y^  will  exceed  Yt  .  To  quantify  this  ,  w>  folio.  Kfroa  [(>  j 
and  utiliz-  th»  Kh  pin  r.-M<  if  r  (til],  [6],  [l])  estimator  of  th°  true 
irvival  function. 


Let  Z,,N  <  •••  <1/  x  denote  th-  order  d  valu  s  of  th-  Z's. 
(1)  -  -  (n) 

Then  th-  non  parametric  maximum  like  lihood  estimator  of  the  survivor 
ship  function  H(t)  -  P(Y  >  tj  for  t  >  0  is  given  by: 


Vi 


H(t)  =  ff  { (n-r )/ (n-r+l) }  , 


vh*  re-  r  •  •  through  tho.,  poaitive  iiitr^<>rr>  for  which  *'(r)  ^  *■ 

nnd  Z,  ^  i;  !IU  uncennomd  observation.  Our  K: i pi m ir — M< •  i  r  adjusted  b 
.•cores  ar-  -riv« -n  in  Tnbl< 


Tabl«  . 

Valuer  of  b.  .  ba 
iJ 

ced  on  th''  Kaplan-M 

1  i  r  .  ntimator. 

(5.  .  ) 

i  . 

a  >  'ii 

Z,  =  Z. 
i  J 

"i  <  ^ 

(1,1) 

1 

0 

-1 

(0.1) 

1 

1 

(H(I1  J)/H(Zi))-l 

(1.0) 

l-  >{h(z1)/k(j'.  )} 

-1 

-1 

(0.0) 

)/il(  .)) 

i-fn(z.)/H(z})} 

fii(  ,)/»( ".)H 

i.'ot’  that .  ;:c  pt  for  tlv  cj:;i  Z.  -  Z,.  the  b.  .  neon  m 

i  J  iJ 


fjiven  by 


bi.i  -  r(yi  >Yjlzi-zJ'Bi-8j’1')-1 


vh'T"  th  eor.iitional  probability  ic  interpreted  ac  if  Y^  nnd  Y 

A 

v:  re  actu;  11,,'  drawn  from  II. 

Tlif  Kaplnn-M-  i  r  ncor>  .  conditional  tent  in  bn,  d  on  0,  and 
is  defined  by  (6),  wh  r  now  tlv  b'r  nn  obtained  row  T-  I >  1  * •  and 

in  l>  etiou  i.  Unucr  Asnuinpt  i  on  (*\) ,  In  t.-nt 


th  a  '  n 


a r*  an  u 
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i  conditionally  dist  ribution- free.  With  the  b  scores  given  by 
T  il-  .  tlv  conditional  moan  and  variance  for  the  nonn  1  approxi¬ 
mation  r  given  by  (7)  and  (8),  respectively. 

t.  A  pti  udo-condit.lonal  tent.  In  Section  3  we  pointed  out  that 
Assumption  (A)  van  not  satisfied  whai  X  =  date  of  transplant  and 
Y  -•  length  of  post , -transplant  survival.  The  hypothesis  of  independ¬ 
ence  of  X  and  Y  in  of  interest  because  a  positive  association 
would  nugg.  st  a  beneficial  trend  in  pout-transplant,  survival  inc< 
initiation  of  th’  .Stanford  Program.  Since  the  procedures  of  Sections 
3  and  •.  are  not  applicabli  in  this  situation,  we  we>v  It  d  to  the 
following  "p.  udo-condit ional"  test. 

When  5.  =  0.  we  know  that  Y.  >  B.,  and  we  also  know  B.  .  We 
1  i  i7  i 

car.  th  ~n  "fill  in"  the  unknown  Y^  value  by  taking  a  random  value 
from  the  distribution  function  defined  by 


(10) 


P(Yi  >  t|Y.  >  Bj)  =  H(t)/H(Bj'  , 


where  H(t)  i.  th  Kaplon-Meicr  er.tim.utor  defined  by  (°).  However, 
the  rand  cur.  ti  tribution  functions  d<  fined  by  (0)  and  (]0)  ar«  not 
continuous  and  there  is  also  ambiguity  in  the  Kn plan-Mi 'i< -r  estimator 
when  Z^j,  the  Largest  Z  vvlu<  ,  corresponds  to  ■  censored  vu3u<  . 

A 

In  that  case,  Kaplan  ar.d  Meier  sugge.  ted  that  for  I  >  Z^j,  ii 

A, 

'.ould  b  i-eg  rad  as  being  b  uw<  n  H(Z^  and  0,  but  not 
r*r.  r  p<  ci  i  d. 

A 

Ir.  ora  r  to  sample  from  II,  and  to  avoid  obtaining  tie  v  lu>  .  ,  we 
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completed  and  adjusted  H  to  H*  (say)  as  follows.  We  formed  a  poly 

A 

onal  curve,  connecting  the  points  (Z^,H(Z^)),  where  r  runs 

through  the  subscripts  corresponding  to  uncensored  values.  We  then 

iefir.ed  H*(t)  to  be  the  polygonal  curve  up  to  Z,  x  ,  and  bevond  ?,, 

\n)  (n 

(11)  H*(t)  =  oxp(Un  fe(n))]tA'.(n))  ,  t  >  Z(n)  . 


This  corr  s ponds  to  fitting  an  exponential  to  the  tail  as  follows. 

For  t  >  ^  y  H*  is  given  by  the  probability  Ii*(t)  =  exp{-t/0}, 

A 

vh  r  *  i.  chosen  so  that  for  t  =  Z yy  cxp{-t/fl)  -  H(Z^),  the 
F.  .plan-!v  i  r  iraute  for  ,t  =  7,y  y  W<>  used  this  adjustment  for  all 

c  s-T.  that  is,  wh  a  corresponded  to  o  censoia'd  value  and  when 

^(n)  corr  o.poah  d  to  an  unc  nsored  value.  Let  Y*  denote  a  random 
value  sampled  from  the  distribution  defined  by 


(12)  P(Y^  >  t)  =  ]I,'(t)/H“(Bi)  , 

where  i  run  • hrough  thorn  subscripts  for  which  =  0. 

Our  pseudo-conditional  test  consists  of  applying  the  usual  per¬ 
mutation  test  based  on  H  to  the  uncensored  "sample"  [(X^,Y^),i 
such  that,  =  1,  (X^,Yp,  i  such  that  =  0] .  If  choosing  the 
Y*  1  by  random  mechanism  is  too  distasteful  to  the  user,  an  nlter- 
n  tiv  ru  -m  .  io  i.  ‘.o  take  Y*  to  br  the  median  of  the  distri¬ 
bution  defined  by  (l  ). 
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6.  Appli cations .  In  this  section  w  give  the  results  ol‘  the 
following  correlation  analyses:  (i)  sex  versus  post-transplant  :  ur- 
vival  time,  (ii)  age  -it  transplant  versus  post-transplant  survival 
time,  (iii)  date  of  transplant  versus  post -transplant  survival  time, 
(dv)  waiting  time  for  n  donor  ve  rsus  post  -transplant  survival  time, 
(v)  date  of  acceptance  versus  waiting  t^rne  for  a  donor,  (vi)  date  of 
•  cot  ptanC'  v  rsur  sex,  (vii  )  date  of  acceptance  versus  age  at  accept¬ 
ance  date. 

The  tests  advocated  in  Sections  3,  and  'j  are  referred  t.o  as 
.  i tuple  adjusted  (JA),  Kapla n-Meier  -adjusted  (KMA),  and  p  eudo- 
conditioual  (PC),  respectively.  Assumpti on  (A)  of  fi.ct.iou  i  is  not 
reasonable  for  projects  (iii)  and  (v)  and  thus  th<  .A  and  KMA  price- 
d  *re  aspect,  in  those  cases;  nevertheless ,  for  comple*  eness ,  w. 

repor  th  "'.d  KMA  result, s.  In  projects  (vi)  u.d  (vii)  then  is 
no  censoring  a •  .'A,  KMA,  and  PC  al.l  reduce  to  the  same  procedure, 

namely  t.l»«  usual  test.  of  i  ndepetdeuce  based  on  Kendall'  !t. 

Our  snr.pl.  sixes,  ire  n  =  (projects  (i)-(iv))  'rnd  n  =  HP 
(project.  -  vii)  .  Thus  n  in  too  large  to  perform  n!  permuta¬ 
tions.  Our  results  are  stated  in  terns  of  the  normal  approximation 
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ar.d  an  approximation  to  the  permutation  ter.t  based  on  a  random  sample 


of  1000  permutations.  In  addition  to  giving  the  one-sided  signi ficanee 


probability  (?  value)  associated  with  S,  we  include  the  value  of 


(13)  !'  =  S/(  £  a*  £  bf  }1/2  , 

i,J  J  i  ?  J  1J 


th«  latter  being  inti  rpn  table  as  a  meacun  of  correlation.  Th 


results  are  summarized  in  Table  3 • 


Table 

3.  Significance  probabilities 

;  for  corr< 

lotion  analyc 

P 

P 

Normal 

(Normal 

(1000  per¬ 

Project 

Procedure 

r» 

deviate 

approx. ) 

mutations  ) 

p 

i 

rA 

-72 

-  .665 

.253 

.?91 

-.080 

i 

KKA 

-62.1 

-.553 

.290 

.  2 

-.067 

i 

PC 

-.446 

.328 

.349 

-.058 

ii 

CA 

-184 

-2.046 

.020 

.0° 

-.205 

ii 

KI-iA 

-521.6 

-  .Ii8 

.016 

.019 

-.  '14 

ii 

PC 

-58  6 

-2.320 

.010 

.008 

iii 

GA 

428 

1.618 

.035 

•  0J3 

.1 33 

iii 

KT, 

.04.1 

I.656 

.049 

.048 

.166 

iii 

FC 

524 

.068 

.019 

.017 

.19fl 

iv 

l'A 

642 

2.729 

.0032 

.003 

•  27  3 

iv 

KK>\ 

657-7 

-.697 

.0035 

<.001 

•  •  1 

iv 

PC 

681 

.701 

.0035 

.001 

.  >0 

V 

ca 

-226 

-.567 

.285 

.277 

-.04  ' 

V 

KMA 

-208.9 

-.457 

.324 

.315 

-.OJo 

V 

PC 

-404 

-.809 

.209 

.187 

-.Ool 

Vi 

IA 

276 

1.078 

.140 

.164 

.O' >8 

vii 

CA 

-338 

-.677 

.249 

.  86 

-.0*  1 

(i)  versus  po:  t-tranr.plant  survival;  H<  r  X  a  ]  if  (lie 

patient  i  mala,  0  if  fi-mal  ,  and  Z  is  the-  uiiniimnn  of  t li<  l.im> 


from  transplant  to  uxith  and  the  time  from  transplant  to  cloning  data  . 
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Eighteen  observations,  corresponding  to  those  patients  who  received 
transplants  1  tore  the  closing  date  and  were  alive  at  the  closing  date 
are  censored.  The  data  are  sparse  for  detecting  sex  differences  since 
only  four  of  the  fifty-two  transplant  patients  are  vomen.  The  data  do 
not  indicate  a  significant  correlation  between  sex  and  post-transplant 

A 

survival.  The  Kaplnn-Moicr  estimator  H(t)  of  P{Y  >  t},  where 
Y  =  post -transplant  survival  time,  is  given  in  Table  4,  In  addition 

A 

to  K(t),  Table  4  contains  an  estimator  of  the  standard  deviation  of 

A 

H(t).  The  estimator,  given  by  Kaplan  and  Meier  [ll],  is 


U~) 


£(fi(t))  =  H(t)(z[(n-r)(n-r+l)]'1}1/ 
r 


where  Z^  <  •••  <  Z,  ^  are  the  ordered  Z's,  and  r  runs  through 
those  posit  iV'  inter  rs  for  which  t  and  ^(r)  an  Uil" 

c or. sored  obC( tv  tion . 

Til  plm-M'  ier  estimator  is  based  on  th<  assumption  that  the 
Y's  (of  ,.*hich  only  34  a r'  observable  her<  )  are  independent  and  id. -r  - 


tically  distribut  'd.  There  is  som-  evidence  (see  project  (ill))  Uni 
in  fact  there  is  a  tr^nd  in  the  Y's  and  so  this  assumption  is  prob- 

A 

ably  not  vnli  .  nevertheless,  we  need  H(t)  in  order  to  apply  the 
KMA  and  PC  procedures.  Furtheimorc,  Table  4  is  ini’onnntivi  in  its 
own  right. 
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T*.  Lie  4.  KnpLan-Meier  estimator  of  the  post-transplant 
survival  distribution. 


t  (day;,) 

H(t) 

SD(H(t)) 

0 

.931 

.019 

1. 

.962 

.027 

3 

.942 

.032 

10 

.923 

.037 

12 

.904 

.041 

15 

.885 

.044 

24 

.865 

.047 

25 

.8iit> 

.050 

39 

.827 

.052 

43 

.808 

.055 

44 

.788 

.057 

.768 

.059 

47 

.749 

.060 

51 

.709 

.063 

54 

.689 

.065 

6o 

.670 

.066 

6i 

.650 

.067 

63 

.629 

.068 

64 

.609 

.068 

65 

.568 

.070 

127 

.547 

.070 

136 

.526 

.071 

147 

.505 

.071 

228 

.482 

.071 

253 

.459 

.071 

31? 

.435 

.071 

322 

.410 

.071 

624 

•  379 

•  078 

730 

.341 

.075 

83  6 

.299 

.077 

1024 

•  >39 

.081 

1850 

.119 

.094 

Apt  at  transplant 

versus  post' 

-transplant  survival:  Hi 

X  =  age  at  transplant,  Z  in  an  in  project  (i).  Procedures  fiA,  KMA , 
ana  PC  all  indent'-  a  small  but  significant  negative  coma  la  ti  on 


i  i  i  mi 
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-  v.."-en  go  at  transplant  and  post-transplant  survival  (see  Table  3). 

(iii)  Pat.'  of  transplant  versus  post-trnnspl  nt  survival:  II  re 
--  date  of  transplant,  Z  is  as  in  project  (i).  Assumption  (A)  is 
..ot  valid,  but  we  upplied  procedures  f<A  and  KMA  to  compare  them  with 
tit  PC  test.  The  analyses  (see  Table  3)  suggest  a  modest  but  signifi¬ 
es-  r.t  positive  correlation  between  post -transplant  survival  time  and 
date  of  transplant,  with  more  recent  candidates  in  the  program  gener¬ 
ally  faring  better  than  earlier  candidates.  One  would  like  to  attri¬ 
bute  this  trend  to  improved  surgical  technique,  and  po -t- surgical 
care,  but  it  may  be  due  to  other  factors  as  well,  .such  as  selection 
of  hardier  piti-ait.s  or  mor  appropriate  candidates  in  recent  years  of 
the  program.  Projects  (v)-(vii)  investigate  whether  there  has  been  a 
change  over  tirm  in  vailing  t  iin  1  rs  a  donor,  sex  of  the  candidate's, 
and  age  of  the  candidates. 


(iv)  iting  time  for  a  donor  versus  post-transplant  surviv.n 

ti.au  :  H  r  -  tiro*  ITom  acceptance  to  transplant,  '  is  as  in 
project  (i).  Proa-  lur<  KMA,  and  PC  all  indicate  ignificani 

positiv'  corr'lation  between  waiting  tiin  and  post -transplant  ur- 
vivul  (so  T  11  3).  On  *  possible  explanation:  Thor  •  *  patients  th.  t 

have  to  «••?.* i"  ■  long  tiro*.*  for  a  donor  (and  survive  to  receive  a  heart.) 
nr*  selected  or  hardiness  by  this  process,  and  might  be  expected  to 
surviv  longei  . 

(v)  Dat«*  o  acceptance  ynimun  waiting  time  or  a  donor;  Hi  n 
X  =  acceptance  cute  and.  Z  is  the  minimum  of  thi  timer'  to  death. 
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transplant,  deselection. 

and  closing  date. 

The  uncensored  values  ] 

are  those  Z's  corresponding  to  patients 

who  received  transplants 

Thirty  of  the  eighty-two 

obsf  rvatlons  are 

censored.  Ther  is  no  e 

dence  (see  Tabl  3)  of  a 

trend  in  waiting  time  for  a  donor.  The 

Ksplnn-Meier  estimator  of  P(Y  >  t),  wher 

e  Y  =  waiting  time  for  1 

donor,  is  giv<  n  in  Table 

9. 

Table  9 .  Knpinn- 

Meier  estimator  o 

f  -uniting  time  for 

a  donor 

•  distribution. 

/N 

A  ^ 

L  (dt|ypJ 

ulli 

SD(H(t)) 

0 

.976 

.017 

i 

.991 

.027 

2 

.912 

.032 

3 

.399 

.037 

7 

.899 

.070 

5 

.8I15 

.071 

7 

.831 

.073 

9 

.803 

.0% 

11 

-776 

.078 

13 

•  799 

.090 

19 

•7J'9 

.091 

16 

•  719 

.093 

17 

.700 

.097 

13 

.689 

.099 

19 

.670 

.096 

20 

.639 

•  097 

22 

.627 

CO 

Lf\ 

c 

• 

27 

.608 

.099 

'25 

•  977 

.060 

26 

.976 

.060 

27 

.530 

.061 

30 

.919 

.061 

31 

.799 

.061 

32 

.783 

.061 

39 

.790 

.061 

36 

•  73 

.061 

" 


I 


A  4  A 


t  (days ) 

H(0 

SD(H(t)) 

.35^ 

.061 

50 

.312 

.001 

56 

.291 

.060 

57 

.270 

.059 

66 

.248 

.059 

70 

922[ 

.097 

77 

.203 

.096 

82 

.158 

.092 

(vi) 

Date  of  acceptance 

versus  s<>x: 

Here  neitlv  r  variable  is 

censored . 

The  independence  to 

st  yields  a 

non-signil'icant,  trend 

nT 

II 

• 

with  the  indication 

of  a  slight 

.  tendency  tow  rd  choosing 

male  candidates  as  the  program 

progresses . 

(vii) 

Dat<  of  acceptance 

versus  age 

at  acceptance  date;  If  iUr  r 

variable  is  censored,  and,  on 

the  basis  of  the  r*.  test,  oiv  i;  led  to 

accept  th 

hypothesis  of  no  association. 

See  Tabl>  3* 

7.  Robustness  contents. 

Th"  DA  and 

KMA  procedures  are  exact 

vh  n  Assumption  (A)  is  satisfied.  The  PC  procedure  is  not,  exact  but 
a  hint  of  i»  robustness  can  be  obtained  by  comparing  it  with  !!A  and 
KMA  for  situations  where  Assumption  (A)  is  valid.  Projects  (i)  and 
(3 i)  furnish  such  comparisons,  and  there*  is  reasonable  agreement  (  e 
Table  •).  For  project  (iv),  when  the  deviation  from  As,  umption  (A) 
is  mild,  2A,  KMA,  and  PC  arc  also  in  coo1  agreement. 

In  project  (iii),  th  deviation  from  Assumption  (A)  i;  strung; 

-1  co  r  lesion  between  X  and  B  exists  (wh'*Ui.  r  or  not  X  u<l 
Y  are  iso*  p'-r.  .  at ) .  Here  th"  significance  probabilities  achieved  iiy 


i 


i 


t 


« 


2Lt 

SA  and  KMA  may  overemphasize  the  signii'i ounce  of  a  positive  correlation 
between  X  and  Y.  The  fact  that  the  P  value  for  PC  is  even 
mailer  than  those  for  SA  and  KMA  may  be  in  part  due  to  the  exponen¬ 
tial  tail  adjustment  (see  (ll),  Section  5)  ol’  the  Kaplun-Meier  esti¬ 
mator.  The  exponential  tail  tends  to  be  optimistic  about  post- 
transplant  survival  prediction — especially  for  those  •>  nsored  values 
corresponding  to  transplants  that  occurred  shortly  befor<  the  closing 
da  te , 

In  view  of  the  importunci  of  project  (iii),  the  following,  mor< 
conserve tiv-  ,  approach  to  post-transplant  survival  prediction  was  also 

A 

used.  After  connecting  the  consecutive  points  (Z^^  ,11(7  ^. ^ ) )  as  in 
Section  f; .  th  line  connecting  the  points  corresponding  to  the  two 
largest  imconsored  values  was  extended  downward  until  it  inter;  ot  d 
th  horizontal  axis.  Ihis  tail  replaced  th<  exponential  tail  given 
by  (ll).  Using  this  straight-line  adjustment,  to  th<  Kapiun-Meii  r 
estimator,  th-  PC  procedure  yielded  a  correlation  of  r  =  .109  with 
an  estimated  one-sided  significance  probability  of  P  =  .112  (bn;  e<! 
on  1000  i  mutations ).  The  corresponding  values,  using  the  exponen¬ 
tial  tail  for  the  Knplnn-Meier  estimator,  were  p  =  .108  and 
?  =  .017 .  Tli'  overall  indication,  ’roin  th'-  PC  procedun  s,  is  that  o 
a  modest,,  but  significant,  positive  correlation  between  transplant 
j  te  .ap/i  post-tr an;, plant  survival. 

We  obtained  another  measure  of  trend  in  port-transplant  survival 
follows.  Table  6,  in  th  form  of  a  2x4  coating- ncy  t  • . i  1  , 
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giver,  for  each  of  the  years  1%8,  1S>69,  1970,  1971>  th<  number  of 
patients  (M)  who  w<  re  accepted  as  candidates  in  that  year  and  went  on 
*.o  r  ceive  a  transplant.  It  also  provid<  s  the  number  (m)  of  those 
who  survived  at  least  one  y.  ur  nl’ter  their  transplant.  (We  omitted 
1-7  c  to  simplify  th<  analysis;  complete  returns  on  the  1972  group 
re  not  available  at  the  March  1,  1973  closing  date.) 


Table  6.  Relationship  between  year  of  acceptance  and 
one-year  post-transplant  survival  rate. 


year  i 

1(1968) 

2(1969) 

3(1970) 

^(1971) 

rai 

2 

6 

2 

6 

M .  -m . 

l  l 

_L 

c 

h 

_L 

o 

11 

6 

13 

Kmploying  a  procedure  due  to  Annitage  [  1  ] ,  we  used  Kendall's  f”. 
to  t.<  ct  for  a  trend  in  the  one-year  survival  rates.  The  statistic 
r  luces  to 

2  =  (M1-m1)(m2+m^m,v)  +  (M  -in,)(m,+m^)  +  (M 
-  rr^f  (M,  -m.)  +  (M^h.)  +  0\-m^)}  -  m 
+  (”)rmj,)}  -  m3(Mu-n>4)  • 

r.i  r  -h  hypothesis  of  no  trend,  the  mean  of  f>  is  zero  and  variance 
of  T  is 


I  1 


h  r  t  =  nnd  T  =  J>1  .  For  the  data  of  Table  6,  S  =  50. 

7  r(f)  =  7539,6,  and  the  normal  deviate  is  .7^  corresponding  to  a 
o:.  -  id-1  P  value  of  .23.  This  analysis,  based  on  much  less  infor¬ 
mation  than  PC  procedures,  agrees  in  direction --but  not  in  degree  of 
sign!  f leaner— with  the  PC  results  that  suggested  a  positive  trend  in 
post. -transplant  survival. 
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